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REFLECTION  MATRIX  FOR  OPTICAL  RESONATORS 
IN  FEL  OSCILLATORS 


I.  INTRODUCTION 

Free  Electron  Lasers  (FEL)  operating  as  oscillators^  ^  require  the 
trapping  of  light  pulses  between  systems  of  mirrois  ( *. esona tors ) .  ^  ^ 
These  pulses  are  repeatedly  amplified  via  synchronous  interaction  with 
electron  pulses  passing  through  the  viggler.  The  simplest  oscillator 
configuration  is  that  of  an  open  resonator  with  two  opposed  identical 
mirrors.  The  radiation  vector  potential  for  this  arrangement  is  expressed 

in  terms  of  the  free  space  eigenmodes  A  (r)  =  A  (r)e  of  the  paraxial 

mn  mn  mn 

wave  equation^,  where  e  is  the  polarization  vector  and 
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mn 


(1) 


0  1 

The  exponent  ^>(r)  2  k[z  +  (xt*y")/2R(z)  ]  in  (1)  contains  the  rapidly 

varying  phase  on  the  wavelength  scale  X  =  2n/k.  The  constant  phase 

2  o 

wavefronts  are  spherical  of  curvature  1/R(z)  -  z/(z^  *  b').  The  spot  size 

221/2  1/2 

of  the  radiation  envelope  is  W(z)  =  v  (1  +  z  /b  )x  ,  where  w  =  (2b/k) 

is  the  waist,  and  the  distance  z  is  measured  from  the  position  of  the 

waist.  The  amplitude  squared  of  the  mode  d’ops  by  1/2  over  a  distance 

equal  to  the  Rayleigh  length  b  (also  known  as  confocal  parameter).  Most  of 

the  radiation  is  confined  within  a  cone  parametrized  by  the  diffraction 

1/2 

angle  9,  =  W/z  =  (X/bn.)  .  The  amplitude  profile  u  (r;V)  contains  the 

a  mn 

transverse  spatial  variation,  equivalent  to  a  small  k,.  perpendicular  to 
the  z-direction.  Higher  modes  correspond  to  a  larger  effective  k  , 
resnonsible  for  the  phase  shift  exp  I i 6  (z)l.  ror  a  given  k  the  mode  is 

completely  defined  by  the  two  independent  parameters  R  and  *  (01  an\ 
combination  of  two  out  of  the  four  quantities  R.  w,  z  and  b). 

The  functions  u  (r;V)  depend  on  the  elected  coordinate  system.  In 
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rectangular  coordinates  (x,y,z)  they  are  given  by 


u  (x , y ; W)  =  a  H  H  f^)  e  2  ,  (2) 

mn  J  mnmfVjnfV/  V 

where  H  ,  H  are  the  Hermite  polynomials,  a  =  <?/U^ )  ^ /^(  ri  2m+n  m!  n!) 
m  n  r  •  mn 

is  the  normalization  factor,  setting  the  total  electromagnetic  flux  carried 
by  the  mode  equal  to  unity,  and  =  (m  +  n  +  1)  tan  \z/b)  is  the  slow 

phase.  In  cylindrical  coordinates  (r,G,z),  where  tan©  -  x/y,  r  = 
(x  +y  )  ,  ur(r,9;W)  take  tne  form 


u~P( r, ©;U) 


'cosp©' 
,s  i  n  p  0, 


(3) 


where  +p  (-p)  signifies  cosine  (sine)  poloidal  dependence,  LP  are  the 
Laguerre  polynomials,  ap  =  (l/2rtW^)^,/‘'  |m!/(m  +  p ) !  ]  and  &p(z)  =  (2m  + 

p  +  1 )  tan  ^ ( z/b) . 

The  electron  beam  is  an  optically  active  medium  that  alters  the 

characteristic  parameters  of  the  radiation  after  each  passage.  During  the 

build-up  period  the  modal  content  and  the  structure  of  the  light  pulses 

inside  the  oscillator  will  change.  A  numerical  method  has  been  developed 

recently  optimizing  the  representation  for  the  amplified  radiation.  In  the 

12  13 

source  dependent  expansion  ^ ’  the  waist  size  and  the  curvature  of  the 
elected  modal  basis  is  tailored  according  to  the  driving  source  term.  That 
minimizes  the  number  of  modes  required  to  describe  the  light  beam.  In 
general,  the  curvature  and  waist  size  of  these  modes  does  not  match  the 
curvature  and  waist  of  the  vacuum  eigenmodes  for  the  resonator.  Therefore 
the  transfer  matrix  for  a  given  mirror  must  be  known  for  arbitrary  incoming 
modes.  This  need  stems  from  computational  as  well  as  physical  reasons. 
The  knowledge  of  the  cavity  reflection  matrix  R,  together  with  the  gain 


matrix  G  th-ough  the  viggler,  is  necessary  in  determining  the  potential  for 
steady  state  opcicttion. 

During  high  power  operation,  grazing  mirror  incidence  may  be  necessary 
to  avoid  exceeding  rhe  dielectric  breakdown  limit  for  the  reflection 
surface.  Also,  in  case  of  a  high  per-pass  gain  with  optical  guiding,  the 
waist  for  the  reflected  radiation  would  be  much  larger  than  the  incoming. 
In  two-mirror  resonators  the  reflected  radiation  could  bit  the  wiggler. 
Therefore,  ring  resonators,  including  three  or  more  mirrors,  must  be 
employed  for  the  deflection  and  recirculation  of  the  light  pulses.  The 
study  of  the  reflection  matrix  must  be  extended  to  include  the  effects  of 
deflecting  the  light  beam,  in  addition  to  finite  mirror  size  and  curvature 
mismatches . 

The  amplitude  profile  of  the  incoming  radiation  will  be  modified  by 
reflection.  A  single  incident  mode  Amn(rj)  will*  in  general,  be  partially 
reflected  into  different  modes  Apq(rf))  where  (m.n)  t  (p,q).  This  is  caused 
by  the  deflection  of  the  light  beam,  the  finite  size  of  the  mirror  and  the 
curvature  mismatches.  Reflection  into  other  modes  will  affect  the 
interaction  between  the  electron  beam  and  the  radiation  in  a  number  of 
ways.  First  the  rms  radius  of  the  light  beam  will  change,  affecting  the 
matching  beam  condition.  Second,  the  light  pulse  will  spread  axially 
because  of  dispersion  among  different  modes,  since  the  phase  velocity 
depends  on  the  modal  number  (m,n).  Third,  different  phase  advances  during 
reflection  among  the  various  modes  may  render  these  modes  out  of  phase 
after  a  number  of  bouncings  off  the  resonator.  For  the  above  reasons  the 
fraction  of  radiation  scattered  into  other  modes  will  contribute  to  the 


losses  in  FEL  oscillators. 


II.  OUTLINE  OF  i HE  METHOD 


The  geometry  of  the  reflection  is  illustrated  in  Fig,  1.  The 
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reflection.  The  origins  are  displaced  from  the  mirror  center  by  1^  and  1 

respectively,  where  1^  is  the  distance  of  the  minimum  waist  w.  for  the 

incoming  radiation  and  1  is  the  distance  of  the  minimum  waist  v  for  the 

o  o 


reflected  mode.  A  third 

coordinate 
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in  the 
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surface  coordinates  in  each  reference  frame.  The  transformations  among  the 
va’-ious  frames  are  defined  by 
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We  are  ’nterested  in  cases  when  the  reflected  radiation  remains 
focused  along  some  direction  zq  making  an  angle  <f>  with  the  incoming  z.. 
Then  the  reflected  vector  potential  will  also  be  expandable  in  free  space 
eigenmodes  Apq(r0)  of  the  paraxial  equation  in  the  new  direction.  The 
mirror  surface  generating  focused  reflection  in  the  desired  direction  can 
not  be  arbitrary  but  must  be  appropriately  defined.  The  angle  of 
deflection  <t>  will  enter  the  equation  defining  the  mirror  surface.  The  other 
surface  parameter,  namely  the  curvature  1/R  ,  is  a  free  parameter.  It 
determines  the  curvature  1/R  for  the  outgoing  modes  given  the  curvature 
1/fL  of  the  incoming  modes.  In  case  of  reflection  by  an  arbitrary  mirror 


4 


surface,  the  scattered  radiation  cannot,  in  general,  be  covered  by  the 

paraxial  modes  that  do  not  form  a  complete  set  in  three  dimensions. 

We  consider  incoming  radiation  of  given  curvature  and  of  arbitrary 

amplitude  profile  A1(r^),  consisting  of  various  modes  (m,n)  with  the  same 

R.(z).  If  both  incident  and  reflected  radiation  are  expanded  into 
1 

eigenmodes , 


A1 ( r . )  =  E  c  A  ( i . ;  , 

'  1 '  u  mn  mn  1 
m ,  n 


A°(r  )  =  E  c°  A  (r  )  , 

°  p,q  Pq  P^  ° 


the  relation  among  the  incident  and  reflected  expansion  coefficients  c1^, 

c°  is  formulated  by 
pq 

C°=  R  c1  ,  (6a) 

or 


O  yn  1 

C  =  >  R  C 

pq  pq  mn 

m ,  n 


(6b) 


where  Rmn  are  the  elements  of  the  reflection  matrix  R. 

pq 

We  examine  the  case  when  the  mirror  dimensions  p  are  much  larger  than 
the  wavelength  X,  X  <<  p  (otherwise  diffraction  rather  than  reflection 
would  prevail).  We  also  assume  that  the  angle  C  subtended  by  the  mirror 
C  =  p/R  ,  where  R^  parametrizes  the  radius  of  curvature,  is  small,  of  the 
order  of  the  diffraction  angle  9^,  C  ~  9^  ~  t.  The  i-th  component  of  the 
reflected  vector  potential  at  distance  |r  -  rQ  |  >>  X  from  the  mirror 

surface  S  is  then  given  by 


5 


In  Eq  .  (7)  n  •  Ar  is  the  obliqueness  factor  vheit  Ar  =  (  rQ-r  ) .' j  r  -  rQ  |  ail(J  11 

is  the  normal  unit  vector  to  the  reflecting  surface.  The  surface  element 

ds  is  given  by  ds  =  &[z  -f(x  ,y  )  1  dx  dy  dz  vhere  z  =  f(>:  ,y  )  is  the 
6  '  1  o  o  -o  ‘  o  Jo  o  o  o  Jo 

3 

surface  equation.  Equation  (7)  is  the  convolution  of  a  source  term  A  (r  ) 
at  the  mirror  surface  with  the  propagator  exp( ik  |  c0~r0  | ) / |  r  - rQ  | ,  i.e..  a 

3 

superposition  of  spherical  waves  originating  at  S.  The  source  term  A'fr^) 
is  specified  from  the  incoming  vector  potential  AJ(ix)  through  the  boundary 
conditions  and  the  coordinate  transformations  (4).  We  will  assume  a 
perfectly  conducting  surface,  where  the  incident  and  reflected  fields  are 
related  by 

As  =  -A*  +  2  (n  •  A*)  n  ,  (8a) 


and  n  is  the  normal  unit  vector  to  the  reflecting  surface.  When  the  solid 
angle  subtended  by  the  mirror  is  small,  the  last  term  in  (8a)  is  very  small 
and  the  boundary  cone  it  ton  becomes, 


s  ,  i 

(V)  =  ’  *  (V) 


(8b) 


Relation  (8b)  corresponds  to  a  phase  shift  by  t  dvr’-'g  Tr  is 

independent  of  the  wave  polarization,  thus  the  subscript  (v)  is  dropped. 

Cross  polarization  effects,  due  to  the  last  term  in  (8a)  are  discussed  in 

Ref.  10.  Most  of  the  computations  will  be  performed  on  the  mirror  surface. 

To  simplify  the  notation  from  now  on,  we  drop  the  bar  (  )  under  the  mirror 

coordinates  r.  Subscripted  quantities  such  as  r^,  rQ,  will  signify  the 

mirror  surface  in  each  reference  frame.  IJnsubscripted  coordinates  will 

denote  the  observation  point  in  the  reflected  radiation  frame  of  reference. 

We  seek  cases  when  the  reflected  radiation  propagates  focused  along 

z  ,  contained  within  a  cross  section  or  dimensions  >:,v  <<  z-z  .  The 

o  o 
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O  o  ? 

expansion  |r-r  )  =  (z-zq)  (1  +  [(x-xq)  +  ( y  -  y  o )  ^  )/2(z-zo)‘‘ }  replaces  the 
full  propagator  inside  (7)  with  the  paraxial  propagator  U  ^  in  that 
d i rec  t ion , 


-ik(z-z  ) 
o 


U  -  ( r ,  r  ) 
-k  o 


(x-xo)2+  (y-y0)2 
2(z-zo) 


It  is  known  that  the  profile  of  a  given  eigenmode  A ( x0 > >’0 * z0 )  at  z0 
is  generated  by  the  propagator  U^(r,rQ)  acting  on  the  mode  A^^fx.y.O)  at 
z=0.  The  inverse  propagator  U  ,  (r,rQ)  therefore  reproduces  Amn(x,y,0)  from 
( Xq . yo , Zq ) •  This  suggests  multiplying  and  dividing  the  integrand  inside 
(7)  by  exp  [  i  *(  r  )  ]  /  [1  -  z  2/bQ2  ]  ,  recasting  (7)  in  the  form, 


(10) 


(ID 


and  the  phase  A(rQ)  is  given,  in  outgoing  coordinates  rQ,  by 


A(  r  )  =  k  z . (r  )  +  zn  + 


(  r  )  +  y  .  (  r  )  xc+v"l 

10  J10  O  0 


2R . ( r  ) 
i  o 


2R  (r  ) 

o  o 


•  (12) 


The  expression  for  A(rQ)  depends  on  the  angle  through  the  transformations 
between  the  incoming  and  the  outgoing  coordinates,  Eqs.  (4).  Expression 
(10)  is  the  approximation  of  the  exact  solution  (7)  to  order 


7 


,  (x-x  )z  »  ( v - v  )‘\/2(z-z  )“  ~  C  .  It  is  valid  provided  that  the  surface  S 
produces  focused  reflection  along  the  desired  direction.  (Uhet vise,  tire 
paraxial  limit  will  fail  to  encompass  all  the  radiation  contained  in  the 
original  expression  (7). 

The  term  exp|ifi(r  )1  is  varying  rapidly,  on  the  scale  of  the 
wavelength  X.  Therefore.  its  convolution  with  the  slowly  varying  source 
term  over  an  arbitrary  surface  will  he  vanishingly  small.  In  general,  this 
corresponds  to  radiation  scattering  where  only  a  small  fraction  of  the 
incoming  radiation  is  reflected  along  the  considered  direction  $.  The 
integral  (10)  will  he  finite  only  when  it  is  possible  to  satisfy  the- 

condition  £>(  r  )  =  constant  over  some  surface  S.  If.  in  addition,  S  is  much 

o 

larger  than  X.  expression  (10)  will  he  finite  only  within  a  narrow  angle  &<J> 
around  <f.  This  guarantees  that  the  reflected  radiation  remains  focused 
along  that  direction.  Therefore,  a  condi ' ion  that  the  exact  reflected 
radiation  (?)  be  fully  coveted  by  the  paraxial  lim.it  (1")  is  that 


£s(  rQ)  =  cons  t  .  ,  (13) 

along  'he  surface  S.  Accordingly.  the  optical  path  is  the  same  along  the- 
various  rays  connecting  an  incoming  wavefront  with  its  mirror  image 
(reflected)  wavefront. 

Requirement  (:;)  defines  (tic-  appropriate  mirror  surface  r  =f  (x  ,y  ; *) 
for  reflection  in  the  elected  direction.  Expressing  all  quantities  inside 
(13)  in  the  mirror  coordinate  frame.  applying  the  transform.. lions  (A)  and 

9 

us i nr  the  scaling  >'  /R  ~  v  /R  ~  z  <<  1,  z  /F.  ~  t  ve  obtain  i tom  (13) 

s  m  ■  s  rr,  s  rr, 


1 


2R  cos 
m 


& 

9 


(Ha) 


where 


X 


(14b) 


1  1  1 

R  =  2R  *'  2R . 
m  o  l 

Equation  (14a)  is  the  analytic  expression  for  a  paraboloid  surface. 
R  parametrizes  the  mirror  curvature,  being  positive  or  negative  for  a 
convex  or  concave  mirror  respectively-  The  surface  is  reflection  symmetric, 
with  \zy.)  ^  and  (zy)  as  the  symmetry  planes;  there  is  no  rotational 
symmetry  around  z  .  Surface  (14a)  can  also  be  approximated,  to  second 

n  o 

oidc:  in  (:<  /R  )*',  (y  /R  )~  by  hyperboloids  or  ellipsoids  defined  by 


-  R  cos 
m 

fl  - 

o  o  .  o 

*i.  <T 

X  COS  ~  -  V 

s  2  •  s 

D  2  2  <t> 

=  R  COS  7T  . 
m  2 

(15a) 

+  R  cos 
rc 

9 

ir  * 

2  2  ♦  2 

X  COS  ~  -  v 

c;  J  *  q 

2  * 

R  cos  -  . 
rc  2 

(15b) 

All  the  surfaces  become  spherical  in  the  limit  of  perpendicular  incidence  e 
=  0,  and  plane  mirrors  when  F.^  -  Relation  (l-+b)  defines  the  curvatuie 

of  the  reflected  modes  from  the  incoming  mode  curvature  and  the  curvature 
of  the  m i r r o  : . 

Switching  Eq .  (12)  into  the  mirror-aligned  coordinates  r  through  Ecs. 
(4),  anc,  using  the  surface  constraints  (14),  it  follows  that 

A(  r  •  const.  -  oj  kp  ,  j  . 

'■'m'  J 

where  p  parametrizes  the  mi  i  rot  size.  A  more  cor.pl  icated  surface  equation 
(higher  than  quadratic  in  ,  y.  z)  is  required  to  improve  the  constancy 
a  higher  order.  Since  kp  >>  1.  the  approximation  A(r^)  const,  is 

satisfactory  for  a  first  order  expansion  of  the  ret  lection  matrix  in  powe: 
of  p/R  .  as  long  as  p/R  <  (kp)  .  In  rase  p  P.  >  (kp)  the  s  i  ■ 

fl!  Hi  —  !'■ 

variation  of  A(x( ,y)  over  S  must  he  included.  That  introduces  or. 
additional  contribution  in  the  reflection  matrix.  known  as  sphetical 
abera  t i on . 
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III.  COMPUTATION  OF  THE  REFLECTION'  MATRIX 
The  reflected  radiation  is  expressed  by 


o ,  , 
t.  (r) 


- 1  4>  ( x  ,  y  ) 

,  ,  ,  .  ov  o  'o'  ,,  ,  V 

dx  dy  o(x  ,y  )  e  U  ,  (r, r  ;  , 

o  Jo  v  o  •  o  -k  o 


(16) 


where  a(x  ,y  )  =  S[x  , y  ,z  (x  , y  )].  Expanding  the  source  a(x  , y  )  in 
o  o  o  o  o  o-'o  1  o-o 

terms  of  u  (>:  ,  y  ), 
mn  o  ■  o 


a(x  ,v  )  -  )  R  u  (x  ,y  ;V  ). 
o  '  o  ^  mn  o  J  o  o 

m ,  n 


(17) 


and  exploiting  the  property  of  the  inverse  propagator  U  ^ ,  the  reflected 


vector  potential  A°(r)  at  z  -  0  becomes 


.  o,  r-.  mn  ,  ..  . 

n  { ,  v  .  0 )  ^  r  F.  u  ( x .  v ; »  ) . 
■  ,-J  mn  -  o 

m ,  n 


(18) 


where  V  (z)  -  w  (1 
o  o 


2..  2.1/2 
?  /  b  ) 
o 


1/2 

(2h  / k )  .  Expression  (18)  is  a 


complete  decomposition  of  the  reflected  radiation  into  paraxial  eigenmodes 
for  incident  radiation  of  arbitrary  profile. 

According  to  the  definition  (6b).  the  R"n  element  of  the  reflection 

pq 

matrix  R  is  obtained  from  the  source  term  c  (:•:  , v  )  inside  (16)  generated 

pq  o  •  o 

try  a  single  incident  eigenmode  A  [ r . ( rQ ) ] .  The  integration  is  performed 
in  the  n-.i  rror  -aligned  coordinates.  taking  advantage  of  the  existing 
symmetries;.  The  coordinates  r^  and  i  ,  defining  the  incoming  and  outgoing 
wave  functions;,  become  explicit  functions  of  >:  ,  y(.  through  the 
t  ratisf ormat ions  (4).  The  surface  equation  (14a)  is  used  to  express  r  in 


terms  of  (:•:  ,y^).  The  mirror  boundary 


to 


2  2  <(> 

>:s  cos  2  +  y< 


is  defined  by  the  intersection  of  the  infinite  surface  (14a)  with  the  plane 

9  9 

z  =  const  =  2  p^cos^f  4>-'2 )/  R  .  After  the  above  manipulations  the 

s  '  m  1 

reflection  matrix  elements  take  the  form 


dx  dy 
s  s 


u  ( x  ,  v  )  u  < x  ,  y  ) 
run  s  s  pq  '  s  s 


2  ,1/2 


z  ( x  ,  v  ) 
O  S  •  s 


Zi  (xs’ys} 


i A( x  ,v  ) 
s  s 


i6x  (x  , y  ) 
pq  s  J  s 


.  V  ) 

■ 

X 

, 

x  sin  s’. 

S  'S 

4>  (, 

s  .  4> 

s  2) 

cos  -  1 

R  Sln  2 
m 

0  1 

*-  \ 

R  "  ; 

m 

where 


u  ( x  ,  y  )  h  u  x  ( x  ,  v  ) ,  y  ,  u  ( x  ,  v  )  =  u  x .  ( x  ,  y  ) ,  v 
mn  s  s  nn[  o  s  ’s  JsJ  pq  s  ■  s  pq  L  l  s  ■  s  ■  sj 

9  2 

Expression  (20)  is  correct  to  order  pVR 

Each  representation  of  R  is  tied  to  the  choice  of  the  basis  function.'; 
u  (r).  In  any  case  R,  as  given  by  (20).  depends  on  four  parameters^ 

R  s  R(4>,  a.  u;  O.  (21) 

if  is  the  reflection  angle  shown  in  Fig.  1.  a  is  the  ratio  of  the  incoming 

to  the  outgoing  spot  size  at  the  mirror,  a  -  V . ( 1 . ) / W  (1  ).  u  =  p/w 

110  0  o 

parametrizes  the  mirror  size  compared  to  the  radiation  spot  size.  E  = 


V  /R  scales  as  the  diffraction  angle  0,  =  V  /]  multiplied  Ivy  the 

o  m  d  o  o 

curvature  mismatch  R  /R  between  the  mirror  and  the  radiation  wavefronts. 

o  m 

The  spot  size  enters  as  a  free  parameter  because  only  the  curvature  1  ■' Rq 
for  the  reflected  modes  is  specified  by  the  mirror  geometry.  Since  many 


combinations  of  Vq  and  1  apply  to  a  given  curvature  according  to  paragraph 

Eq .  (1),  an  additional  selection  rule  for  V  is  needed.  Note  that  V  does 

not  have  to  match  W . .  This  is  obvious  in  cases  when  the  mirror  size  p  is 

smaller  that  V. .  Each  value  of  V  defines  a  complete  set  of  modes  for  the 
i  o 

reflected  radiation  and  an  equivalent  representation  for  R. 

Parameters  <f>,  a,  and  u  can  be  arbitrary.  In  most  cases  of  interest, 
however,  E  is  small,  E  <<  1,  of  the  same  order  as  the  diffraction  angle  9^. 
The  analytic  computation  of  the  matrix  elements  is  carried  out  by  expanding 
the  integral  (20)  in  powers  of  E, 

R  =  R(0)  +  E  R(l)  +  E2  R( 2 ) .  (22) 

The  first  order  expansion  is  performed  in  Ref.  10.  In  this  paper  we 
review  some  of  the  general  properties  of  R  and  focus  on  the  reflection 
of  the  lowest  mode  Uqq. 


IV.  LIMITING  CASES 

When  the  mirror  radius  tends  to  infinity  (1/R  -»  0),  or  in  cases 

m 

of  vertical  incidence  on  the  mirror  ( 4>  =  0),  the  higher  order 

corrections  in  the  reflection  matrix  R  disappear, 

R  =  R(0)  (23) 

in  both  representations.  The  nondiagonal  elements  in  R  stem  from  the 
finite  mirror  size  only.  If,  in  addition,  the  mirror  size  is  very 
large ,  u  >>  1,  it  is  appropriate  to  take  =  V .  as  best  representation 
for  the  reflected  radiation.  The  a  =  1  limit  yields 


R 


mn 

pq 


,mn 

6 

pq 


(24) 


Thus,  in  case  of  large  curved  mirror  and  vertical  incidence,  or  large 
plane  mirror  and  arbitrary  incidence,  the  reflection  matrix  is  the 


12 


identity  matrix. 

The  case  a  =  1  is  of  special  interest  for  arbittary  angle  of 
deflection  4>  and  mirror  curvature  1/R,  as  it  will  be  explained  in  the 
next  section.  For  finite  mirror  size  p  >  W  ,  (  u  >  1),  there  exists 
l  zeroth  order  nondiagonal  terms  inside  R ( 0 ) .  Since  R(0)  is  independent 

of  the  angle  of  deflection  <f>,  the  mirror  size  yields  the  dominant 
contribution  to  the  reflection  into  modes  different  than  the  incoming. 
The  effects  of  the  deflection  of  the  light  beam  enter  to  first  order  in 
E,  R(l),  or  higher. 

As  the  mirror  size  becomes  very  large  and  the  limits  of  integration 
are  extended  to  infinity  the  orthogonality  among  the  various  modes  u 
becomes  effective.  The  off-diagonal  terms  in  R(0)  become  comparable  to 
the  first  order  terms  R(l)  roughly  when  l/u"  ~  E  ~  0^  .  At  the  limit  U 
®  all  the  nondiagonal  elements  of  R  are  reduced  to  order  E  or  higher. 

Rmn  =  £  Rmn(1)  ,  0(E2).  m  *  p,  n  *  q.  (25a) 

pq  pq 

and  the  only  matrix  elements  of  zeroth  order  in  E  are  the  diagonal 

Rmn  =  Rmn  Q  „  0(£,2).  (25b) 

mn  mr. 

in  both  Hermite  and  Laguerre  represenat i ons .  The  lowest  correction  in 

2 

the  diagonal  elements  is  of  second  order  V~ ,  while  the  first  order 
contribution  disappears.  This  is  consistent  with  flux  conservation 
during  reflection  in  case  of  a  large  mirror. 

The  superposition  principle  can  be  used  to  describe  reflection 
from  more  complex  mirror  surfaces.  In  case  of  a  mirror  with  a  hole  the 
surface  integral  (7)  over  S  is  expressed  as  J<.  -  where  is 

defined  by  the  mirror  surface  including  the  hole  surface,  and  Sp  is  the 
surface  of  the  hole  only.  The  total  reflection  matrix  R  is  given  by  R  -- 


13 


R(S.)  -  R(S0),  the  diffetence  in  the  reflection  matrices  associated  vith 
mirrors  and  S2  respectively.  The  transmission  matrix  T  through  a 
screen  with  an  aperture  of  area  S  is  given  by  T  =  -  R,  R  being  the 
reflection  matrix  for  a  mirror  matching  the  aperture  S.  The 
transmission  matrix  for  radiation  diffracted  behind  a  finite  size  mirror 
is  given  by  T'  =  1  -  elfl  R  where  1  is  the  identity  matrix. 


V.  REFLECTION  OF  THE  LOWEST  ORDER  MODE 

The  computation  of  all  the  truncated  integrals  for  finite  mirror 

suiface  is  nontrivial.  Most  appl i _a t ions  however  involve  the  (0,0) 

lowest  order  mode  as  the  dominant  mode  in  both  incoming  and  reflected 

radiation.  The  strategy  here  is  to  compute  the  element  R^  of  the 

reflection  matrix  first.  Then  the  waist  for  the  reflected  modes  W  can 

0 

be  selected  so  that  it  maximizes  R^J.  The  optimum  representation 
cond i t ion 


9R 


00 

00 


da 


=  0, 


(26) 


puts  the  maximum  amount  of  the  reflected  radiation  in  the  lowest  mode  (a 
different  mode  and  matrix  element  may  be  chosen,  if  desired).  It  is 
pointed  out  that  (26)  does  not  improve  the  properties  of  the  reflected 
radiation.  It  enables  one  to  choose  the  best  representation  in  terms  of 
minimizing  the  coefficients  of  the  undesired  modes  for  the  scattered 
radiation.  Once  W  is  fixed  by  (26)  then  the  exact  location  and  size  of 
the  waist(s)  for  the  reflected  modes  is  determined  by  solving  the  system  of 
equations 


U 

o 


2,  1/2 


o 


l 

R 


1 

o 


o  o 


(27) 
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The  element  is  identical  in  both  representations  since  the  lowest 
order  mode  is  the  same  in  rectangular  and  cylindrical  coordinates. 
Performing  the  integration  (20)  yields  to  first  order  in  f 


Note  that  the  first  order  term  vanishes  and  the  lowest  correction,  is  of 

2 

second  order  in  £,  .  The  exact  dependence  on  the  mirror  size  p  is 
parametrized  by  u  =  p/VQ,  while  a  =  V . /Wq  parametrizes  the  ratio  of  the 
incoming  and  scattered  radiation  spot  sizes  at  the  mirror.  The  optimization 

00  n  o  n  /  o  o 

condition  3R0()(0)/3a  =  0  yields,  a"  =  l  *  exp| -( 1  +  )p“  ]  (  a  +( 2p  +1 ) a  -1  ]  . 

In  case  the  mirror  cross  section  is  much  larger  than  the  spot  size  of  the 
incoming  mode,  p  >>  1,  a  -*  1  and  the  reflected  spot  size  at  the  mirror 
matches  the  incoming,  V °  =  V. . 

Large  mirror  size  is  desired  to  maximize  the  total  reflection 


coefficient.  The  reflection  coefficient  hR  is  given  by  =  P  /P.  where 

'  i  1 2  i  i  2 

incoming  flux  is  =  jc  \  =  L  |c  )**  and  the  outgoing  flux  is  given  by 


and  the  outgoing  flux  is  given  by 


E  E  IR 

mn  pq 


mn 

pq 


2 


(29) 


In  Fig.  2  we  plot  r^  for  the  lowest  order  incoming  mode  as  a  function  of 
U'  =  cos  (<(>/?)  p/V  =  cos(<f>/2)  u  ■  v'  parametrizes  the  size  of  the  mirror 
surface  projection  into  the  plane  perpendicular  to  the  incoming 

radiation  direction.  The  incoming  radiation  has  a  wavelength  X  -  Ip 

_4  _  2  o 

(10  cm),  waist  v.  =  10  cm  at  distance  1.  =  2.5x10  cm  from  the  mirror 
i  i 

and  radius  of  curvature  (at  the  mirror)  =  2.5x10  cm.  The  mirror  has 

a  radius  of  curvature  -  2.5x10  cm,  yielding  reflected  modes  of  R^ 


)5 


3 

2.5x10  (again  1  and  vq  depend  on  the  choice  of  Vq).  In  Fig.  3  we  plot 

the  reflection  coefficients  of  the  lowest  order  mode  (0,0)  into  the 

Pd 

first  25  modes  (p,q)  with  p  <  q  <  5,  as  a  function  of  u'.  The 

deflection  angle  is  90°  and  the  ratio  of  the  spot  sizes  is  1. 

Increasing  mirror  size  maximizes  the  diagonal  element  and  minimizes 
scattering  into  other  modes. 

In  Fig.  4  we  fix  the  mirror  size  u'  =  2  and  the  angle  4>  =  90°  and 
vary  the  spot  size  ratio  a.  The  best  representation,  maximizing 
and  minimizing  is  obtained  at  a  =  1.  However,  for  small  mirror  p '  = 

0.66,  the  maximum  for  R^  occurs  at  a  =  0.70  (see  Fig.  5).  Radiation 
reflected  off  mirrors  smaller  than  the  incoming  spot  size  is  best 

described  by  outgoing  modes  of  reduced  spot  size  <  V . .  Also  note 
from  Fig.  5b  that  for  small  mirror  size  the  total  power  reflected  into 
the  first  25  modes  never  exceeds  80°„  of  the  incoming  flux;  even  with 
many  more  modes  remains  less  than  1.  In  Fig.  6  the  reflection 
coefficients  R^  are  plotted  as  functions  of  the  angle  of  deflection 
for  fixed  a  =  1,  u'  =  2.  It  is  seen  that,  for  sufficiently  large 

reflecting  surface  and  good  choice  of  the  spot  size  ,  the  reflection 
matrix  is  not  very  sensitive  to  4>.  Comparison  of  Figs.  2  and  3  with  the 
rest  of  the  plots  shows  that  the  relative  mirror  size  to  the  radiation 
spot  size  is  the  most  important  parameter  to  determine  the  reflection 
into  other  than  the  incoming  modes. 
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